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ABSTRACT 

Analytical expressions for covariances of weak lensing statistics related to the aper- 
ture mass Map are derived for realistic survey geometries such as SNAP for a range 
of smoothing angles and redshift bins. We incorporate the contributions to the noise 
due to the intrinsic ellipticity distribution and the effects of finite size of the cata- 
logue. Extending previous results to the most general case where the overlap of source 
populations is included in a complete analysis of error estimates, we study how vari- 
ous angular scales in various redshifts are correlated and how the estimation scatter 
changes with survey parameters. Dependence on cosmological parameters and source 
redshift distributions are studied in detail. Numerical simulations are used to test the 
validity of various ingredients to our calculations. Correlation coefficients are defined 
in a way that makes them practically independent of cosmology. They can provide 
important tools to cross-correlate one or more different surveys, as well as various 
redshift bins within the same survey or various angular scales from same or differ- 
ent surveys. Dependence of these coefficients on various models of underlying mass 
correlation hierarchy is also studied. Generalisations of these coefficients at the level 
of three-point statistics have the potential to probe the complete shape dependence 
of the underlying bi-spectrum of the matter distribution. A complete error analysis 
incorporating all sources of errors suggest encouraging results for studies using future 
space based weak lensing surveys such as SNAP. 

Key words: Cosmology: theory - gravitational lensing - large-scale structure of 
Universe - Methods: analytical, statistical, numerical 



1 INTRODUCTION 

Detection of weak lensing signals by observational teams (e.g., Bacon, Refregier & Ellis 2000, Hoekstra et al. 2002, Van 
Waerbeke et al. 2000, and Van Waerbeke et al. 2002) has led to a new avenue not only in constraining the background dynamics 
of the universe but in probing the nature of dark matter and dark energy as well. To do so one compares observations with 
theoretical results obtained from simulations or analytical methods. 

Often numerical simulations are employed which use ray-tracing techniques as well as line of sight integration of cosmic 
shear (e.g., Schneider & Weiss, 1988, Jarosszn'ski et al., 1990, Wambsganns, Cen & Ostriker, 1998, Van Waerbeke, Bernardeau 
& Mellier, 1999, and Jain, Seljak & White, 2000, Couchman, Barber & Thomas 1999) to test analytical calculations of cosmic 
shear. Of course simulations are being constantly updated to match upcoming surveys with broader sky coverage. 

Analytical techniques too have made progress in understanding of shear induced by line-of-sight density inhomogeneities. 
Typically cosmological perturbation theory is employed at large angular scales (e.g., Villumsen, 1996, Stebbins, 1996, 
Bernardeau et al., 1997, Jain & Seljak, 1997, Kaiser, 1998, Van Waerbeke, Bernardeau & Mellier, 1999, and Schneider et 
al., 1998) while techniques based on the hierarchical ansatz provide a good match to simulation results at smaller angular 
scales (e.g.. Fry 1984, Schaeffer 1984, Bernardeau & Schaeffer 1992, Szapudi & Szalay 1993, 1997, Munshi, Melott & Coles 
1999, 1999a, 1999b). Ingredients for such calculations include Peacock & Dodds (1996)'s prescription (see Peacock & Smith 
(2000) for a more recent fit) for the evolution of the power spectrum or equivalently the two-point correlation function. Recent 
studies report an excellent agreement between analytical results and numerical simulations of weak lensing effects (Valageas 
2000a & b; Munshi & Jain 2000 & 2001; Munshi 2000; Bernardeau & Valageas 2000; Valageas, Barber & Munshi 2004; Barber, 
Munshi & Valageas 2004; Munshi, Valageas & Barber 2004). One the other hand, one may also use halo models for specific 
purposes (e.g. Takada & Jain 2003). 
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In realistic survey scenarios, weak lensing signals are accompanied by various noises due to the intrinsic ellipticity distri- 
bution of galaxies, shot-noise due to the discreet nature of the source galaxies and finite volume effects due to finite survey 
size. A detailed analysis of these sources of noises was carried out by Munshi & Coles (2003) following Schneider et al. (1998). 
Recent studies by Valageas, Munshi & Barber (2004) have used such techniques to model errors in future surveys such as 
SNAP. It was also possible to study cross- correlation among various surveys using such a technique. Recent analysis by Mun- 
shi & Valageas (2004) showed that cross-correlation studies between two different surveys will have the potential to detect 
systematics associated with weak lensing measurements. In Munshi & Valageas (2004) we concentrated in dividing the source 
population into two non-overlapping subsamples. Generalising these studies here we are also able to probe estimation error 
associated with more than two different smoothing angular scales for same or different source redshift distributions. Besides 
the generalisation of our studies to three-point statistics means that wc arc able to probe the whole bi-spectrum and to use 
it to constrain theories of primordial non-Gaussianity or non-Gaussianity induced by gravitational clustering. 

This paper is organised as follows: in section 2, a very brief introduction to our notations is provided. The window 
functions used in the rest of the paper arc introduced along with a simple model of correlations hierarchy. In section 3, 
we take into account these realistic sources of noise to study the scatter associated with various estimators. Specific survey 
characteristics based on SNAP class experiments are assumed for such an analysis. In section 4 we check the validity of our 
results against numerical simulations for a range of redshift distributions and smoothing angular scales. We also investigate 
the dependence of two-point and three-point correlations on large-scale structures. Finally section 5 is left for a discussion of 
our results. Appendix A provides details of the analytical results used in the main sections. 



2 NOTATIONS AND FORMALISM 

Any weak lensing effect X smoothed over some angular radius 6s and averaged over the redshift distribution of sources n{zs) 

can be written in terms of the fluctuations of the density field as: 

X = Jm Ux0) dx wix) Six,m, with wix) = ^ J^"^^^ dz. n{z,) ^ ^^^^^1^^ (1 + (1) 

Here the redshift z corresponds to the radial distance x ^^d T) is the angular distance, i9 is the angular direction on the sky, 
5{x,'D'&) is the matter density contrast and hereafter we normalize the mean redshift distribution of the sources (e.g. galax;ies) 
to unity: J dzs n{zs) = 1. We note amax the depth of the survey (i.e. n{zs) = for Zs > ^max). Here and in the following we 
use the Born approximation which is well-suited to wcak-lensing studies: the fluctuations of the gravitational potential are 
computed along the unperturbed trajectory of the photon (Kaiser 1992). We also neglect the discrete effects due to the flnite 
number of galaxies. They can be obtained by taking into account the discrete nature of the distribution n{zs). This gives 
corrections of order 1/N to higher-order moments of weak-lensing observables, where A'^ is the number of galaxies within the 
circular field of interest. In practice N is much larger than unity (for a circular window of radius 1 arcmin we expect N > 100 
for the SNAP mission) therefore in this paper wc shall work with eq.(l). The angular flltcr Ux depends on the weak lensing 
observable one considers. For instance, the fllter associated with the aperture-mass Map is (Schneider 1996): 



where G is a top-hat with obvious notations. The angular radius 9s gives the angular scale probed by these smoothed 
observables. As described in Munshi et al. (2004), the cumulants of X can be written in real space as: 



or equivalently we can write in Fourier space: 



/•Xmax P f P f ^ \ 

<^'')= = y^ l[dxiw{xi) J l[dk^Wx{k±jV^es) (n^''""'M Wki)...<5(M)c. 



(4) 



We note (..) the average over different realizations of the density field, ^j, is the real-space p— point correlation function of the 
density field ^p(xi, ..,Xp) = ((5(xi)..(5(xp))c, is the component of k parallel to the line-of-sight, kx is the two-dimensional 
vector formed by the components of k perpendicular to the line-of-sight and WxQix'DOs) is the Fourier transform of the 
window Ux'- 



WxO<i±T>es) = j dSUx{&)e''^-^^, whence VFMap(kx©6l,) = 



24J4(fcj,Pge) 

{k±vesr ■ ^ ' 



The real-space expression (3) is well-suited to models which give an analytic expression for the correlations ^p, like the 
minimal tree-model (Valageas 2000b; Bernardeau & Valageas 2002; Barber et al. 2004) while the Fourier-space expression (4) 
is convenient for models which give a simple expression for the correlations (<5(ki)..<5(kp))c, like the stellar model (Valageas 
et al. 2004; Barber et al. 2004). In this article we shall use the latter stellar model defined by: 
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(5(ki)..5(kp))e = ^fo(ki + ...+kp) ^[]p(fc,), with & = !, (6) 

where Sd is the Dirac distribution and P{k) is the 3-d power-spectrum of the density fluctuations. The coefficients 5*3, 5*4, .. 
are closely related (and approximately equal) to the skewness, kurtosis, .., of the density field. Then, substituting eq.(6) into 
eq.(4) yields (Munshi et al. 2004): 

(XOc = Sp j^"^^^ dx j d# Ux0) Ix{x,^y-\ with 7x(x,^) = \j^ ^-fr^ Wx{\^^Pies) e"*-^^: (7) 

Here we introduced the power per logarithmic wavenumber A^(fc) = 47rA;^P(A;). Of course, these expressions generalize in a 
straightforward fashion for many-point cumulants. For instance, for two-point cumulants we obtain within this model: 

(xfx|)e = Sp+q Z'""" dx wlwi fdS\-^ mi^) hixJ)"-' HxJy + U20) h{xJr /2(x,^)'-'| ■ (8) 

Jo J ~^ 9 P + 1 J 

These quantities describe the cross-correlations between two surveys or subsamples, which are refered to by the subscripts 
"i=l,2". The source rcdshift distributions ni{zs) and the smoothing angles 9si can be different, as well as the directions on 
the sky. However, in this paper we shall restrict ourselves to correlations between angular cells which are centered on the same 
direction. Then there is no need to add an angular shift to the expressions (2) and (5) of the filter functions. On the other 
hand, note that the signal would decrease for nonzero angular separations. In this article we shall consider both two-point 
and three-point cumulants. The latter contain valuable information about the bispectrum and can be used as a template for 
non-Gaussianity studies. 

The weak lensing effects associated with different redshift bins or angular scales are correlated since their lines of sight 
probe the same density fluctuations at low z (where they are largest). In order to measure these cross-correlations we define 

as: 

{XfX|)c . ,. , {XlX2)c /„N 

111 particular m = . „ ' , /. , i /2 ■ (9) 



The quantities Vpq correspond to the two-point cumulants (Xf X|)c normalized in such a way that most of the dependence on 

cosmology and gravitational dynamics cancels out. Thus, if the two subsamples are highly correlated we have rpq ~ 1 while 
rpq ~ if they are almost uncorrolatcd. In a similar fashion we can generalize Vpq to three-point objects rpqs as: 

_ (X^XIXI), . (XlX2X3)c , . 

^^"^ ^j^P+9+f^p/(p+9+'')^j5fP+<i+<'^9/(p+9+'')^XP+<i+<'^V(p+9+'') ' iiiparacuiar rm ^j^3^i/3^j^3^i/3^j^3^i/3 " ^^^> 

These objects provide valuable information about the bi-spcctrum as they directly probe three different angular scales over 
three different redshift distributions. They can also be used to study cross-correlations among three different surveys some of 
which may have overlapping source distributions. However note that as a result of subdividing the same source samples into 
various redshift bins the shot-noise will increase. In addition, being a true three-point object it will be more affected by the 
finite size of the survey volume at larger angular scales. However a good sky coverage for future surveys such as SNAP will 
mean these objects can be measured over reasonable angular length scales. It is possible to investigate even higher numbers 
of points (samples) but we shall not consider this here. 



3 FUTURE SURVEYS 

3.1 A specific case study: Wide SNAP survey 

In order to illustrate the measurements one can expect from future surveys, we consider here the case of the Wide SNAP 
survey. We use the characteristics of the SNAP mission as given in Refregier et al. (2004). The redshift distribution of galaxies 
and the total surface density rig of usable galaxies are: 

n{zs) oc Zg e'^-^'^""'' with zo = 1.13, Zmax = 3 and rig = 100 arcmin"^. (11) 

The shear variance due to intrinsic ellipticities and measurement errors is cr» = (|e,p)'/^ = 0.31 while the survey covers an 
area A = 300 deg^. Therefore, we take for the number N of galaxies within a circular field of radius Os and for the number 
Nc of cells of radius 9s : 

TV = ngTT^? ~ 314 f ^ and iVc = --^ =2.7x10^ f—l^V-A^) (12) 

^ VlOOarcmin-V Vlarcmm/ (26*^)^ V^OOdegV Vlarcmm/ ^ ^ 

For some purposes we shall also consider the two subsamples which can be obtained from the Wide SNAP survey by dividing 
galaxies into two rcdshift bins: Zs > z, (which we refer to as "Widc>") and z < z, ("Widc<"). We choose z, = 1.23, which 
corresponds roughly to the separation provided by the SNAP filters and which splits the Wide SNAP survey into two samples 
with the same number of galaxies (hence rig = 50 arcmin"^). 

For the background cosmology we consider a LCDM model with Qm = 0.3, = 0.7, Hq = 70 km/s/Mpc and as = 0.88. 
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Figure 1. Correlation coefficients rn (left panel), r2i (middle panel) and ri2 (right panel) are plotted as a function of the smoothing 
angle 0s2- AH cases correspond to the full SNAP survey. The three different curves correspond to the smoothing angle 6si = 5', 10' and 
20', whereas 9a2 runs from 1' to 20'. Errorbars associated with each curve show the l-cr scatter in the measured values, due to intrinsic 
ellipticities and cosmic variance. 
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Figure 2. Correlation coefficients rn, r2i and ri2 are plotted for the wide SNAP survey with redshift binning. Solid curves denote the 
complete source sample whereas other line-styles correspond to redshift binning. The dashed curve ">" shows the case where Mapi at 
fixed 6gi is measured from high-2 galaxies (z^i > 1.23) while Map2 at 8^2 which runs from 1' to 20' is measured from low-2 galaxies 
(232 < 1.23). The dotted curve "<" corresponds to Zsi < 1.23 and Zs2 > 1.23. 



3.2 Two-point correlations for Map 

In this article we shall restrict ourselves to the statistics obtained from the aperture-mass Map, defined from the filters Um^j, 
and WAfap given in eqs.(2), (5). We first show in Fig. 1 the lowest-order two-point correlation coefficients Vpq obtained at a 
fixed smoothing angle dsi while 9a2 runs from 1' up to 20', for the full wide SNAP survey. The three curves in each panel 
correspond to 6si = 5', 10' and 20'. The error-bars take into account the cosmic variance and the intrinsic ellipticities of 
galaxies, see Appendix and Munshi & Valageas (2004). 

The correlation coefficients Vpq show a peak where both scales coincide and are fully correlated (r^, = 1 as Mapi and 
Map2 are actually identical at this point). We find from rn that the correlation remains strong ( > 0.4) up to an angular scale 
ds2 twice larger than the other smoothing scale 6si. We must point out that this depends on the smoothing window and that 
top-hat windows (for the smoothed convergence its, or the shear components 7is) generate more extended correlation profiles. 
This is related to the fact that the aperture- mass Map is a very localized filter in Fourier space (Schneider 1996). However, 
we can see that the third-order correlations r2i and ri2 decrease much more slowly (especially r2i). It is interesting to note 
that the shapes of r2i and ri2 are not identical and that r2i exhibits in particular a very extended tail at larger angles. This 
can be understood from eq.(4). The Dirac factor 5D(k-i + ^2) contained in the mean (5(ki)5(k2))c involved in rn implies 
fci = k2 so that rn peaks at 9si = 9s2 where both W{kiT>9si) and W{k2T>0s2) can be simultaneously maximized. Then, rn 
decreases for different angular scales where one cannot simultaneously maximize W{ki) and W{k2). On the other hand, for 
the third-order correlation r2i the Dirac factor is fo(ki -I- k'l -I- k2) with a weight W {ki'DOsi)W {k'i'D9si)W {k2T>9s2) ■ Then, 
at large 9s2 it is still possible to maximize all three weights W by choosing k2 ~ 1/T>9s2, ki ^ k'l ^ 1/T)9si ^ ^2, with 
k'l ~ — ki and k2 = — ki — k'j. This is not possible for ri2 where the weight is now W{k\T>9s\)W{k2'D9s2)W{k'2T)9s2) and 



Covariance of Weak Lensing Observables 5 




I I I I I I I I I I I I I I I I 
5 10 15 20 

B 




5' 



0.8 
0.6 




0.4 - ^ 

0.2 ^ 

--_ ACDM J 

^ ri I I I I I I I I I I I I I I I I r 

5 10 15 20 



Figure 3. Correlation coefficients rn, r2i and ri2 are plotted for various cosmological or survey parameters. The solid line is our baseline 
LCDM model, the "f2ni" curve shows a 20% increase of Om (from flui = 0.3 up to Qm = 0.36), the "erg" curve is a 20% increase of erg 
(from (Tg = 0.88 up to erg = 1.06) and the "2:9" curve is a 20% decrease of the characteristic redshift zq of the survey (from zq = 1.13 
down to ZQ = 0.9), eq.(ll). The distance between these latter curves and the baseline model has been multiplied by a factor 20 since the 
dependence on cosmological or survey parameters actually nearly cancels out. 



from two small wavenumbers k2 and k'2 of order 1/T>6a2 it is impossible to construct a wavenumber ki = — k2 — k2 which is 
much larger (fci is of order 1/T>9s2 <^ 1/T>6si)- A similar reasoning also explains the steep decline of the correlation coefficient 
r2i at low angular scales 9s2 - Thus, the correlation coefficients Vpq directly probe the features of the matter density field (note 
that the previous arguments only depend on the statistical homogeneity of the density field). We can see from Fig. 1 that the 
error-bars expected from the wide SNAP survey are quite small for rn. They are larger for the third-order statistics r2i and 
ri2 but one should still be able to obtain a good measure of the shape of these correlations. On the other hand, the measure 
of the coefficients rpq can provide a way to check the importance of noise and systematics. 

Next, we show in Fig. 2 the effect of the redshift binning of the data. We divide the galaxies into two different bins: 
Zsi > 1.23 and Zs2 < 1.23 (dashed curve ">") or Zsi < 1.23 and Zs2 > 1.23 (dotted curve "<"). Of course, the behaviour of 
the various correlation coefficients rpq remains similar to the case of only one redshift bin displayed in Fig. 1. However, at the 
peak 9si — 9s2 we now have rpq < 1 since even for the same smoothing angle the quantities Mapi and A/ap2 are not identical 
as their lines of sight extend to different redshifts. Of course the error-bars are somewhat larger when we apply such a redshift 
binning of the data but it is still possible to measure these coefficients to a reasonable accuracy up to third-order statistics 
r2i and ri2. This probes the evolution with redshift of gravitational clustering, although the latter is best measured with 
one-point statistics (most of the time dependence of the amplitude of density fluctuations and of the cosmological distances 
associated with various source redshifts is factorized out of the coefficients rpq). 

Finally, we also display in Fig. 3 the results obtained with a 20% increase of f2m (from Jim ~ 0.3 up to Qm = 0.36), or 
a 20% increase of the normalization erg of the density power-spectrum (from erg = 0.88 up to erg — 1.06), or a 20% decrease 
of the characteristic redshift zq (eq.(ll)) of the survey (from zq — 1.13 down to zo — 0.9). The distance between these 
latter curves and the baseline model (solid line) has been multiplied by a factor 20. Thus we see that the dependence on 
cosmological or survey parameters almost completely cancels out of the reduced correlation coefficients Vpq. Therefore, these 
quantities provide a direct handle on the noise or on the detailed properties of the large-scale density field (e.g. the validity 
of a hierarchical model like (6) for many-body correlation functions). On the other hand, cosmological parameters may be 
measured from one-point cumulants. Thus, the use of one-point, two-point and three-point cumulants provides a convenient 
separation between various features of the problem: cosmology, large-scale structures and noise. 



3.3 Three-point correlations for Map 

We show in Fig. 4 the three-point correlation coefficient rm defined in eq.(lO), for the full wide SNAP survey. This three- 
point object is a direct probe of the underlying bi-spectrum of the matter density field. The error-bars are computed from 
the expressions derived in the Appendix where the optimised estimators Hpqs are used (this helps to reduce the scatter, see 
Valageas et al. 2004). The statistic rm is a symmetric function of its arguments when the three source redshift distributions 
are the same as in Fig. 4 where the full source sample is used. Then, when two smoothing angles are equal it coincides with 
the two-point statistics r2i or ri2. As in the case of two-point correlators the error-bars are dominated by the finite size of 
the survey and are higher for larger smoothing angular scales. However from the figure we can see that SNAP will be able to 
probe rm reasonably well even when a large smoothing angle 6s2 = 20' or — 20' is used. Redshift binning can be useful in 
studying the redshift evolution of rm but this reduces the number of sources thereby increasing the shot noise which starts 
playing a dominant role at smaller smoothing angular scales. The noise level at given smoothing angular scales is not only 
determined by the signal-to-noise ratios pi , p2 and ps of the variance of the aperture- mass at each scale but also by the overlap 
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of various smoothing patches as described by pi2, p23 and pis, see eq.(A4) and the Appendix. It is possible to generahse our 
results to cases where the patches are not all centered on the same line of sight but we shall not investigate this point here. 

In Fig. 4 we plot rm for three combinations of {6si,9s2) while runs from 1' to 20'. The curves show a single peak 
reached when 9a3 is equal to the smallest angular scale among {9si,9s2)- This may be understood as follows. Let 9si < 9s2, 
then as O^i increases from 0' the value of rm grows as the surface which is common to the three patches increases as 
^53. More to the point, the relative importance of the common surface increases: relative to the three quantities Mapj it 
goes Eis {(^ss/^si)^, (^33/^32)^, 1}. Next, as 0s3 gets larger than ^si the common surface sticks to (it no longer grows) 
while the relative common surface actually decreases as {1, {9si/9s2)^ , (^si/^ss)^}- Therefore, one can expect a steady decline 
beyond 9si. This is indeed what we observe in Fig. 4. Of course, matters are actually a bit more intricate as the three-point 
correlation is not "^D-correlated" in real space so that the correlation between different patches is not strictly proportional to 
their geometrical overlap and the filters Um^,, are not simple top-hats. However, these simple arguments explain the overall 
trends seen in Fig. 4. 

As expected, for a fixed angle 9si (here 5') the peak of rm at 6^3 — 9ai is larger when 9a2 is closer to ^ai (here 10' as 
compared with 20'). Note that rm — 1 for 9si — 9s2 = 9s3 where the three quantities Mapi, Map2 and Map3 are identical 
whence fully correlated (for identical source distributions). We can see that the correlation of third-order statistics remains 
significant over a large range of angular scales and shows a slow decline at high 9s3 and a steeper falloff at low 9s3, in a fashion 
similar to r2i shown in Fig. 1, for the same reasons. 



4 COMPARISON WITH NUMERICAL SIMULATIONS 

We have shown in the previous section that the two-point and three-point correlators rp, and Vpqs associated with the aperture- 
mass Map could be measured in future surveys like SNAP and allow a clean separation between cosmological parameters and 
the properties of the density field (the behaviour of its many-body correlations). They can also be used to monitor the noise of 
the survey. In this section, we address the use of these correlation coefficients from the point of view of large-scale structures. 
Thus, we first compare our predictions from the simple hierarchical model (6) with the results of numerical simulations. 
Secondly, we investigate the dependence of these correlators on the detailed properties of the density field. 

The numerical method for the computation of the lensing statistics is based on the original formalism of Couchman, 
Barber & Thomas (1999) which concentrates on computation of three-dimensional shear matrices along lines of sight through 
the linked simulation volumes. This technique has been further developed by Barber (2002). This procedure which has been 
described in various publications in detail (e.g. Barber 2002) has been applied to a LCDM cosmological simulation created by 
the Hydra Consortium ^ using the 'Hydra' A''-body hydrodynamics code (Couchman, Thomas & Pearce, 1995). Its cosmological 
parameters are: Qrn = 0.3, Qa ~ 0.7, F = 0.25, Ho = 93 km/s/Mpc and as = 1.22. The computations of correlation functions 
were done by using concentric annuli which are divided in equal radial bins and angular bins in real space directly from the 
rectangular grid. Various combinations of bin width in radial and angular directions were considered while computing the 
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Figure 5. The two-point correlation coefficients rii,r2i and ri2 are plotted for various angular scales. In each plot the curves from left 
to right correspond to 6si = 5', 10' and 20', while 9s2 runs from 1' to 20'. Source redshift is fixed at Zsi = Zs2 = 1. Solid points are 
measurements from simulation data while the solid curve is the theoretical prediction from the simple stellar model (6). Error-bars are 
computed from 10 identical realizations. Note that the LCDM cosmology used for the simulations is slightly different from the one used 
for the SNAP survey. 
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Figure 6. The two-point correlation coefficients rii,r2i and ri2 are plotted for the angular scale 9si = 5', 10' and 20' while 6s2 runs 
from 1' to 20'. The redshift slices considered are z^l = 1 and Zg2 = 2. The solid curve is the analytical model (6) while solid points are 
measurements from simulation data. 

correlations to check the stability of our scheme. Several levels of dilution were tested by increasing the grid size of the parent 
rectangular. For simplicity, we consider Dirac delta distributions of sources. That is, for three-point correlators all sources are 
located at the three redshifts Zai, Zs2 and Zss which can be different. 

4.1 Two-point correlators 

We again study the two-point objects such as Vpq by keeping one angular scale fixed, say 6si, while the other smoothing angle 
6s2 varies from 1' to 20' (unless p = q and Zsi = Zs2, rpq is not a symmetric function of its arguments 9si and 632)- This range 
is determined by the finite resolution of the simulations: smaller scales are affected by the rectangular grid (i.e. resolution) 
while large angular scales are affected by the finite simulation box. Note that effects due to the finite size of the box will be 
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Figure 7. The three-point correlation coefficient rm as a function of smootiiing angle ds'j, for a fixed pair of (Bs\ , 6s2)- The pair {0si , 6s2) 
for each curve is indicated in the plot. The three source redshifts are equal: Zsi = 2s2 = ^sS = 1- The solid curve is the analytical model 
(6) while solid points with error-bars are measurements from simulation data. 



more prominent at smaller redshifts while effects due to the grid play a more important role at larger redshifts. In most of 
our studies we analyse cross-correlations and auto-correlations of various redshift slices against the slice Zs = 1 which is least 
affected by such artifacts and of great observational interest. 

We display the case of identical source redshifts Zsi — Zs2 = 1 in Fig. 5 and the case of different source redshifts in Fig. 6, 
with Zsi = 1 and Zs2 = 2. As in Figs. 1, 2, the two-point correlators always show a peak where both angular scales coincide and 
complete overlap is achieved to generate a high correlation coefficient rpq. Again, in cases where the same redshift slices are 
used for the two smoothing angles the coefficients rpq reach a value of unity, whereas for different redshift slices the maximum 
is less than unity as only partial overlap of the lines of sight is achieved. We also recover the slow decrease at large angles 
ds2 of r2i. We can see that the simple hierarchical model (6) agrees reasonably well with the numerical simulations, although 
there is a small discrepancy at large angles 9a2 for r2i. However, the error-bars in this domain are large. In agreement with 
the nearly total independence of these correlation coefficients onto cosmology, described in Fig. 3, we have checked that an 
OCDM simulation (again from the Hydra Consortium) yields the same results. The error-bars shown are computed from a 
large number of realizations riroai ~ 10. As expected they are larger for higher-order statistics (r2i and ri2 as compared with 
rii). However, note that all realizations we have used are actually constructed from the same N-body density field and are 
not completely independent so there could be a systematic deviation (although expected to be small) in our results which can 
only be studied by using completely different simulations where density fields probed by various lines of sights are completely 
independent. In this article we only consider angular windows with zero separation. Then the scatter is mostly determined 
by the largest of the smoothing angles which defines how many independent patches there are for correlation studies. 

The effects of small box size not only appear as an increased scatter for the estimation of high-order correlation functions 
it can also bias the mean. Detailed schemes exist in case of galaxy clustering to correct or estimate such bias from galaxy 
surveys or simulated catalogs. For large smoothing angles with 6s2 = 20' it is expected that such effects will be visible in the 
computed values of ri2 estimated from our simulated maps as there are only few independent patches within the survey and 
ri2 puts more weight on 6s2 as compared with r2i. Clearly such effects will be more pronounced at larger Oai and at lower 
redshifts where gravitational clustering has effectively generated correlated patches on the sky. We have made no attempt to 
correct for such a deviation. However, note that the compensated filter associated with Map, which is very localized in Fourier 
space, makes these problems less acute than for simple tophat filters. 



4.2 Three-point correlators 

We show in Fig. 7 the three-point correlation coefficient rm as a function of smoothing angle for three fixed pairs of 
{631,632), with equal redshifts Zsi — Zs2 = Z33 — 1. Of course, we recover a behaviour similar to that of Fig. 4 obtained for 
the SNAP survey. Moreover, we can check that the simple analytical model (6) shows a good agreement with the numerical 
simulations. We present in the four panels of Fig. 8 the dependence on the source redshift Zs3 of the correlator rm, with two 
different redshifts Zsi — 1 and Za2 ~ 2. The dependence on Zas is very weak over the range Zsi < Zss < Zs2 while rm decreases 
for larger Zs3, as expected. The agreement of the simple stellar model with numerical simulations is rather good for all cases. 
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Figure 8. The three-point correlator rm is plotted as a function of smoothing angle 8s3 for the fixed pair {dsi,6s2) = (5', 10'). The 
source redshifts are 2:31 = 1 and Zs2 = 2 while we have 233 = 1 (upper left panel), ZsS = 1.5 (upper right panel), 2^3 = 2 (lower left 
panel) and 2^3 = 3 (lower right panel). The solid curve is the analytical model (6) while solid points with error- bars are measurements 
from simulation data. 



4.3 Constraints on the density field 

We have seen in Fig. 3 tliat the correlation coefficients rpq (or Vpqs) are almost independent of cosmological parameters. 
Therefore, they can be used to monitor the noise of the survey or to constrain the properties of the matter density field (e.g. 
the angular behaviour of the many-body correlation functions). We have checked in Figs. 5-8 that the simple stellar model 
(6) (which is a peculiar case of the more general class of hierarchical models) agrees rather well with numerical simulations. 
We now show that this agreement is not due to the independence of the correlators Vpqs on the underlying density field. That 
is, although most of the dependence on cosmology has been factorized out of the ratios (9) and (10) as well as the overall 
amplitude of the density fluctuations (both properties are mostly probed by the one-point statistics), some information on 
the detailed angular behaviour of the density correlations remains in the correlation coefficients rpqs- 

To investigate this point we compare in Figs. 9-10 the predictions obtained from different analytical models. The solid 
line (A) is the stellar model (6) which we have used so far. Here we must note that this model is not fully defined by eq.(6) 
because the parameters Sp vary with scale and time. More precisely, the skewness parameter 5*3 is obtained by interpolating 
between the quasi-linear theory prediction and the HEPT approximation (Scoccimarro et al. 1998), see Munshi et al. (2004). 
Then, 5^3 depends on the slope of the linear power-spectrum Phik) and on the amplitude of the non- linear power A^(fc) at the 
wavenumber k = A/{T>6s). However, for the two-point correlations we have two angular scales 6si, Os2 which can be different. 
Then, we extended the stellar model (6) to these cases by using: 

model (A): Sp,q = S;i^^+'^ S^,/^l+''\ (13) 

where we note Sp^q the coefficient Sp+q which appears in a cross-product of the form (XfX|)c, and the coefficient 5p+p,o 
(resp. So,p+q) is associated with only one angular scale 831 (resp. ds2) where there is no ambiguity for the typical wavenumber 
k. We have checked in sections 4.1, 4.2 that this simple model (A) agrees quite well with numerical simulations. Next, in order 
to evaluate the sensitivity of our predictions onto the approximation (13) we define the variant (B) by: 

model (B): Sp^q = mm{Sp+q.o; So,p+q)- (14) 

Thus the two models (A) and (B) coincide for 9si — 8s2 and both can be described as stellar models of the form (6). The 
freedom between models (A) and (B), due to the scale dependence of the parameters Sp, expresses the fact that the model 
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Figure 9. The two-point correlation coefficients rpq as in Fig. 5, for Qs\ = 5', 10' and 20' and Zsi = Zs2 = 1- Tlie different line-styles 
represent different analytical models: the basic stellar model (A), the variant (B) (see main text) and the Dirac-product model (D). 



(6) is a simple phenomenological prescription but not a fully consistent model (this would require a constant parameter Sp). 
Finally, we also define a third model (D) where many-body density correlations reduce to products of Dirac factors in real 
space: 

model (D): {(5(xi) . . . (5(xp))c = (5d(x2 — xi) . . . (5d(xp — xi). (15) 

We take the normalization coefficient in eq.(15) to be constant so that it cancels out of the correlation coefficients rpq. For 
p = 2 the model (D) merely corresponds to a white-noise power-spectrum, while for all p > 2 it corresponds to the stellar 
model with a white-noise power-spectrum. This model is not expected to be a good approximation of the matter density field. 
Its only purpose is to check whether the detailed behaviour of large-scale structures can be seen in the correlation coefficients 
rpq. It gives the contribution to the coefficients rpq provided by sheer geometrical factors, since we have for identical source 
redshift distributions: 

model (D): (XfX|)e oc J UiiSfUi^)" if ni{zs) = n2{zs), (16) 

where the normalization factor which only depends on (p + q) cancels out of the coefficients rpq. The models (A), (B) and (D) 
extend in a straightforward manner to three-point correlators like rm. Although the models (A), (B) and (D) yield different 
results despite their similarity, we can note that the predictions obtained from model (A) for the third-order correlators we 
study here are actually common to all models which obey {(5(ki)5(k2)(5(k3))c cx SdO^i +k2 +\i.3)[P{ki)P{k2)F{\ii, U.2) +sym.] 
where _F'(ki,k2) is any function of the angle between ki and k2. Indeed, the angular dependence contained in F(ki,k2) 
is integrated out in the same manner for one-point, two-point and three-point cumulants so that it cancels out of rpqs. Of 
course, the behaviour of i^(ki, k2) would show up if we consider different smoothing cells which are not centered on the same 
line-of-sight. 

We show in Figs. 9-10 these three models (A) (solid line), (B) (dashed line) and (D) (dot-dashed line). We can see that the 
differences are quite apparent and are larger than the error-bars obtained from numerical simulations in Figs. 5-8. Therefore, 
although the correlation coefficients rpq and rpqa are almost insensitive to cosmology they show a significant dependence on 
the properties of the many-body density correlations (note that models (A) and (B) are quite close). This separation between 
the background cosmology (which can be probed by one-point statistics) and the large-scale structures themselves should be 
quite useful to discriminate between both classes of phenomena and to probe gravitational clustering. In particular, we can 
note that the differences between various models seen in Figs. 9-10 are sufficiently large to be detected by future surveys like 
SNAP, see the error-bars in Figs. 1-4. Finally, it appears that the simple stellar model (A) defined from eqs.(6),(13) works 
best and provides a good description of the actual density field for these purposes, as checked in sections 4.1, 4.2. 



5 DISCUSSION 

We have extended the formalism developed in Munshi & Coles (2002) to study the cross-correlations among various angular 
scales covering a range of redshift bins. We define the correlation coefficients rpqs in such a way as to minimize the dependence 
on cosmological parameters and on the redshift distribution of the sources. Since the aperture-mass window is very narrow 
in Fourier space they provide a useful tool to study gravitational clustering, through both the power spectrum and the 
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Figure 10. The three-point correlation coefficient rm as in Fig. 7, at x^i = Zg2 = •2^s3 = 1 for three pairs of smoothing angles {0si, Qa2)- 
The different line-styles represent different analytical models: the basic stellar model (A), the variant (B) (see main text) and the 
Dirac-product model (D). 



bi-spectrum. In particular, higher order coefficients r^qs where p -I- g -I- s > 2 carry useful information about the level of non- 
Gaussianity in the projected density field. Generalizing our previous studies from two-point statistics to three-point statistics 
we have used rm as a diagnostics of the bi-spectrum. Indeed, although one-point or two-point cumulants or correlators such as 
r2i also contain useful information about the bi-spectrum, rm allows one to probe the bi-spectrum in more details. These tools 
can be used not only to probe gravitational instability scenarios but also to put constraints on primordial non-Gaussianities 
or higher-dimensional theories of gravitation. We plan to study such issues in more detail in future works. On the other hand, 
these estimators can also be used to monitor hidden systematics. 

We have presented a detailed analysis of these correlation coefficients for the wide SNAP survey, as an illustration for 
future weak-lensing surveys. We have computed both the expected signal and its error-bars, taking into account the intrinsic 
ellipticity of galaxies and the cosmic variance. The cumulant estimators Hpqs are built from the moment estimators Mpqs 
so as to lower their scatter. Note that we take into account the non-Gaussianities of the density field in the computation 
of the error-bars. The scatter associated with SNAP class experiments is extremely small for second-order statistics such as 
the variance or the correlation rn but it increases very fast with the order of the statistics. Nevertheless, these surveys are 
still able to probe third-order statistics with a reasonable accuracy. This should allow one to get meaningful constraints on 
large-scale structures or on the observational noise. 

Next, we have compared a simple stellar model with the results obtained from numerical simulations. We had shown 
in previous papers that this simple analytical model agrees reasonably well with simulations for one-point cumulants (e.g., 
Munshi et al. 2004). We have found in this work that it also gives good predictions for two-point and three-point statistics 
over the entire range of redshifts and smoothing scales of interest. Therefore, this simple model provides a useful tool to 
investigate in great details weak-lensing effects measured through smoothed statistics like the aperture-mass. 

Finally, we have checked that the agreement between this simple model and numerical simulations is not due to a weak 
sensitivity of these correlators onto the structure of the underlying density field. Indeed, introducing two other models for the 
detailed behaviour of the matter density correlations we have found that the correlation coefficients rpqs show a significant 
dependence on the properties of large-scale structures, which can be detected in future surveys like SNAP. Therefore, such 
observations can yield constraints on the structure of the matter density field. Note however that at order three the predictions 
obtained from the simple stellar model (6) are shared by a wider class of models (e.g. all usual hierarchical models). In order to 
discriminate these various models one would need to go up to fourth-order statistics (where one can draw several hierarchical 
diagrams for instance) but error-bars become too large to give good constraints. 

In this paper we have only considered the cases where the smoothing patches are all lined up along the same radial 
direction, although they can have different redshift distributions or smoothing radii. It is trivial to generalise our results 
to cases where various patches are separated by a given separation angle. However with increasing separation the signal to 
noise will degrade which may prevent practical applications. We have only used the compensated filter associated with the 
aperture-mass which can directly be constructed from shear maps. Its presents the advantages of probing a narrow range of 
wavenumbers and of giving a non-zero signal for third-order statistics. Nevertheless our formalism is completely general and 
we can replace this compensated filter by any window (like a top-hat for the convergence or a modified top-hat with angular 
dependence for shear components). However, in such cases the number of patches which can carry completely independent 
information will decrease (because of the contribution of long wavelengths) which will increase observational error-bars. 
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APPENDIX A: 3PT CORRELATORS 

As in Munshi & Valageas (2004), the three-point cumulants {M'^^^M'^^2^ip3) c can be estimated from the estimators Mpqs 
defined by: 

o {Ni-.N-y-.N^,} 

Mpqs — (JVl) {N2) {Nz) ^ ■ . . C^ipCtip C^jiCtji ■ . . f^jg^tjg C^dftii ■ • • ^h^tU- (Aij 

{(il,...,ip);(il,...,jp);(il,...,!p)} 

Here, A^i, N2 and N3 are the number of galaxies within the patches i^O'ii , 7r^^2 a^nd ttS^s for the surveys (1, 2, 3), while Nij and 
are the number of common galaxies to two or all three samples, and we noted {N)p = N{N— 1)..(N —p+1) = N\/{N—p)\. 
In eq.(Al) we only sum over combinations of distinct galaxies and we assumed that Nj ^ 1. Thus, for Mm we sum over 
N1N2N3 - N12N3 - N13N2 - N23N1 + 2iVi23 ^ N1N2N3 terms. We note Qj = Qm^^^j) the value of the filter Qm^^W at 
the location i&j of the galaxy j and e^j its tangential ellipticity. The filter QjWap is given by: 

Ml^, J), 
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The mean and the dispersion a of Mp^s are: 



(A3) 



The scatter a{Mpqs) only involves the variance ai = (|e*|^) of the galax;y intrinsic ellipticities (even if the latter are not 
Gaussian). The latter enters the dispersion through the combinations pj and pij defined by: 



Pj = 



5jV ,.(Mfp,.)e 
3a2 



and pij 



5NiNj(M^piM^pj)c 



(A4) 



The quantity pj is the signal-to-noise ratio of the weak-lensing variance over the random intrinsic ellipticity. It also measures 
the relative importance of the cosmic variance as compared with the galaxy intrinsic ellipticities to the scatter of relevant 
estimators (intrinsic ellipticities can be neglected li p^ 1). 

The estimators Mpqs defined in eq.(Al) correspond to a single direction onto the sky, around which are centered the three 
smoothing windows of radii dgj- In practice, we can average over such directions on the sky with no overlap. This yields 
the estimators Mpqa defined by: 



Mpqs = -r^^M'p'^l, whence {Mpqs) = {Mpqs) = {Ml^^Ml^^Ml^^) and a{Mpqs) = 



n=l 

f{") :„ 



a{Mpqs) 



(A5) 



where Alpqi is the estimator Mpqs for the cell n and wo assumed that these cells are sufficiently well separated so as to be 
uncorrelated. The estimators Mpqs and A4pqs provide a measure of the moments of weak lensing observables. However, as 
shown in Valageas et al. (2004b), it is better to first consider cumulant-inspired estimators Hpqs and Tipqs- For the case of the 
quantity (MapiMap2Map3)c the relevant estimator Hm is: 



1 . -V 

= Mill - A^oilMioo - A^lOlMoio - AllloMool, = ]y- X! '^"l ' "^^^^^^ = (iWaplMap2Map3)c. (A6) 

n=l 

Then the dispersion a^{Hiii) reads: 

CT^(i/lll) = {Ml.Ml^Mls), + (Ma%Mfp3)c(Mfpi), 



1+1 

pi 



+{MliMl^)4Mls)c 



1 + — 

P3 



+ (MaplMap2Map3)c(Map2Map3)c2 



+ (MaVMa^p3)e(Ma^p2) 
1 



1+1 

P2 



1 + 



P23 



+ (MaplM42Map3)c{MaplMap3)c2 



1 + 



Pl3 



+ (MaplMap2Mfp3)e(MaplMap2)c2 



1 + 

Pl2 



+ 2(MaplMfp2)c(MaplMfp3)e + 2(MfpiMap2)c(Map2Mfp3)c 



+2(MfpiMap3)c(M42Map3)c + 3(MaplMap2Map3)c + (M^pl)c(Map2)c(Mfp3)c 
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while the scatter of the estimator Mm can be obtained from: 

Cr^(Mlll) = a^{Hlll) + 2(AffpiMap2A'/ap3)c(Map2Map3)c + 2{Mapl Ma'p2 Afap3 ) c (Mapl Map3 ) c 



+2(MaplMap2Mfp3),(MaplMap2)c + (Mfpi),(Map2Map3)c 



1 + 



Pl 



+ (Map2)c(MaplMap3)c 



P2 



+ (Mfp3)c{MaplMap2)c 



1 + 



P3 



+ (Mapl Map2 ) c (Mapl Map3 ) c (Map2 Map3 ) c 



2 2 2 

6+ + + — 

Pl2 Pl3 P23 



(A8) 



In order to obtain the dispersion of the correlation coefficient rm we also need the scatter of the cross-product cr^(f/iii; -^300) = 
(ii"iiiJ?30o) — (i?iii)(-ff30o) which reads: 



(7li^lli;H300) = (AfaplMap2Map3)c + (MapiAfap2)c(AfaplAfap3)c3 



! + ■ 



Pl3 



+ (MfpiMap3)c(MaplMap2)c3 
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(A9) 
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while (T^(Miii; M300) writes; 



a (Miii;M3oo) 



a^(_fflli;/f300) + (M4l)c(Afap2Map3)c + (M;'piAfap2)c(MaplMap3)c + (M^pi Map3) c ( Mapl Map2 ) 



1 + 



Pl 



+ (M^pi ) e (Mapl Map2) c (Mapl Map3> c 



6+ ■ 



pl2 Pl3 



(AlO) 



We can check that as for one-point and two-point cunmlants, the estimators Hpqs have a lower scatter than their counterparts 
Mpqs and their dispersion only depends on the combination (1 + 1/p). 

Finally, the three-point correlation coefficient rm can be estimated fi-om the estimator 7?.iii: 

'^^111 



^1/3^1/3^1/3 

''•300 ''•030 ''•003 



with (7?.iii) ~ nil 



(All) 



and: 

a'(7^lll) 



2 I cr^ (7^111 ) 1 cr^ (7^300 ) , 1 cr^('Wo3o) , 1 f^l'Woos) _ 2 cr^(Wiii; T^aoo) _ 2 cr^(?^iii;7^o3o) 

''"'I {Wiii)2 "^9 (7Y3oo>2 "^9 (7^o30>2 "^9 (7^oo3>2 3 (7Yiii>(?^300> 3 {Hm) {Hozo) 

2a^(Hiii;Hoo3) _|_ 2 ct^(H3oo; 7^030) _|_ 2 ct^(H3oo; 7^003) _|_ 2 (T^(Ho3o; "^003) I 
3 (Tim) (7^003) 9 (7^300) (7^030) 9 (7^300) (7^003) 9 (7io3o) (7^003) j 



(A12) 



Here we assumed that the scatter of 7?-iii is small and can be linearized. The various terms in eq.(A12) can be obtained from 
the expressions derived above or from Munshi & Valageas (2004). 



